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ABSTRACT 

We report on a new phenomenon of 'alignment' of supergranulation cells in the polar regions of the Sun. 
Recent high-resolution datasets obtained by the Solar Optical Telescope onboard the Hinode satellite 
enabled us to investigate supergranular structures in high-latitude regions of the Sun. We have carried 
out a local helioseismology time-distance analysis of the data, and detected acoustic travel-time variations 
due to the supergranular flows. The supergranulation cells in both the north and south polar regions 
show systematic alignment patterns in the north-south direction. The south-pole datasets obtained in 
a month-long Hinode campaign indicate that the supergranulation alignment property may be quite 
common in the polar regions. We also discuss the latitudinal dependence of the supergranulation cell 
sizes; the data show that the east- west cell size decreases towards higher latitudes. 

Subject headings: convection — Sun: helioseismology 
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the higher-latitude regions, the structure of the cells looks 
more complicated and less organized. 

We take advantage of the Hinode high spatial resolu- 
tion, and study the near-polar regions during the periods 
of the highest inclination of the solar axis to the ecliptic. 
Our finding from the Hinode data analysis is that the su- 
pergranulation cells in the high-latitude regions seem to 
align predominantly in the north-south direction. This 
cannot be explained by the giant-cell scale structures in 
the numerical simulations. 

2. OBSERVATIONS 

The Solar Optical Telesco pe fSOT^lTsuneta et al.l l2008ft 
onboard the Hinode satellite (|Kosugi et al.ll2007t) has made 
several observing runs for local helioseismology. In par- 
ticular, the polar region of the southern hemisphere was 
observed on 2009 March 7, and the polar region in the 
northern hemisphere was observed on 2009 September 25. 
Since the inclination of the solar rotation axis attains its 
maximum (about 7 degrees) in March and September ev- 
ery year, these are the best periods for observing the high- 
latitude regions. For comparison and testing, we also ob- 
tained datasets of a region close to the east limb on 2009 
September 27, and a region around the disk center on 2009 
December 11. All the observations were 16 hours long. 

In March of 2010, SOT carried out nine helioseismology 
observing runs for the south polar region approximately 
once every three days. The observing period of each run 
was 12-16 hours, aside from the observation on March 16, 
which was 6.7 hours long. 

The data acquired with a cadence of 60 seconds were 
Ca II H line intensity maps as well as Fe I 557.6 nm in- 
tensity maps. Since the field of view of the Fe-line maps 
obtained by the Narrowband Filter Imager (NFI) of the 
SOT is wider than that of the Ca-line maps obtained by 
the Broadband Filter Imager (BFI), we set up the tele- 
scope pointing in such a way that the solar limb is in the 
NFI field of view. The NFI data were then used for track- 
ing the limb position and for determining the coordinates 
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1. INTRODUCTION 

Supergranules (jHartl Il954t iLeighton et al.l Il962f ) are 
thought to be convection cells of the Sun with the size of 
20-30 Mm and the lifetime of ~ 1 day. Although more 
than half a century has passed since they were found, 
their physical property is still puzzling and controversial. 
Sim on fe heighten! (|1964l ) suggested that a convective in- 
stability associated with the He ionization zone m ight pro- 
duce the supergranular cells. IStein et all (|2007l ) argued 
that there is nothing special about the supergranulation 
scale except that this is the scale favor ed by magne ti c field s 
interacting w ith c onvective flows. iGizon et al.l ((2003), 
ISchoul (|2003h . and iGreen fe Kosovichev! (|2007t ) suggested 
that supergranul ation has p r opert ies o f traveling convec - 
tion waves, while iRast et al.l <|2004 and lLisle etahl (|2004 
discussed that the supergranulation patterns may be orga- 
niz ed by giant cells. For recen t discussions, see the review 
bv lRieutord fe Rinconl (|2010D . 

Hydrodynamic processes, including supergranulation, in 
the polar regions of the Sun most likely play an impor- 
tant role in the so lar dynamics and m agnetic activity cy- 
cles. For example. iDikpati et al.l (|2010D discussed that the 
speed of the meridional flow in the polar regions could de- 
termine the length of the activity cycle. The processes 
in the polar regions, however, have not been fully inves- 
tigated from the observational point of view, mainly be- 
cause observing the high-latitude regions from the ecliptic 
plane is difficult due to severe foreshortening. Numeri- 
cal simulations of the turbulent convection on the Sun 
have suggested some kind of latitudinal dependence of 
the supergranulation patterns, although convection cells 
of the supergranulation scale have not been clearly repro- 
duced in the global simulations of the co nvection zone, 
excep t in a shallow spherical shell study (|De Rosa et al.l 
2002J). The cells resolved in the global simulations of solar 
convection are mostly of a giant-c ell scale. According to 
the recent simulati on results (e.g.. iBrun fe Toomrdl2u0l 
iMiesch et al.l |2006[) . the convective cells in the lower lat- 
itude regions may align parallel to the rotation axis. In 
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on the Sun in the BFI images. For the helioseismology 
analysis we used the Ca n H data series. 

3. DATA ANALYSIS 

We first aligned the images and determined the helio- 
graphic coordinates of each image. Then, we chose sev- 
eral points on the Sun within the SOT field of view, as 
the central points for the Postel's projection (see below). 
We tracked th e central points with the local differential 
rotation rate (jSnodgrass! [T984) . and remapped each im- 
age into the heliographic coordinates by using the Postel's 
projection. Using the tracked and re mapped images, we 
carried out the time-distan ce analysis (jDuvall et al.lll997t 
iKosovichev fc Duvalllll997ft to measure variations of acous- 
tic travel times and detect subsurface flow structures. 

The SOT observations of the regions close to the limb 
were carried out with fixed telescope pointings. Although 
the i mage stabilization system, the correlation tracker 
(CT;[Sh lmizu et al.l [2008').' was turned on, the field of view 
was gradually drifting toward the disk center mainly due 
to the limb darkening. Therefore, we aligned the images 
by ensuring that the limb stays in the same position in the 
field of view, and then determined the solar coordinates of 
each image. 

When we observe regions near the solar limb, fore- 
shortening is significant, and appropriate remapping is re- 
quired for the data analysis. We us ed the Poste l 's az- 
imuthal equidista nt projection (e.g., iBogart et al.l Il995t 
IZaatri et al.| [20081 . The projected map is what we would 
observe by looking directly from above the region. For 
the dataset taken at the disk center, however, we did not 
project the images, because the distortion effect is negligi- 
ble there. Also, since for the disk center dataset the SOT 
observed the region tracking its center with the rotation 
rate of the Sun, we did not apply any further tracking. 

For the time-distance analysis, cross-correlation func- 
tions of the surface wavefield of the 5-min oscillations 
were calculated from the tracked and projected datasets. 
We measured the acoustic (phase) travel times for several 
distances by fitting a Gabor-type wavelet t o the cross- 
correlation functions (jKosovichev fc Duvalll 1 1997T) . For 
each distance A, we averaged the oscillation signals over an 
annulus of the radius A around a target point, and then 
cross-correlated the averaged signals with the oscillation 
signal observed at the target point. We obtained outward 
(from the target to the annulus) and inward (from the an- 
nulus to the target) travel times at each point in the field 
of view for a set of travel distances A. The difference of 
the outward and inward travel times indicates the diverg- 
ing/converging flow around the target point in the region 
of wave propagation below the surface. Note that for im- 
proving the sig nal-to-noise ratio w e applied phase-speed 
filters following |DmalLetaL| ([1997J). 

4. RESULTS 

4.1. Travel-time maps 

The outward-inward travel-time difference maps of the 
north polar region are illustrated in Figure [TJ Panels (a), 
(b), and (c) show the travel times for the distances, A, of 
8.4, 14.4, and 21.1 Mm, (or 0.7°, 1.2°, and 1.7°) respec- 
tively. The black regions, where the outward travel times 



are shorter than the inward travel times, indicate diverg- 
ing flows, while the white regions correspond to converging 
flows. The dark cells with white boundaries are supergran- 
ular cells. In the maps for the largest annuli (panel c), the 
supergranular patterns are rather faint. Judging from the 
penetrating depths of the acoustic rays, this indicates that 
the supergranules are not much deeper than, say, 5 Mm, 
which is consistent with other hel ioseismology results (e.g., 
IWoodardl [20071 ; ISekii et al.ll2007l ). Hereafter, we use only 
the travel times for A = 14.4 Mm, because these travel- 
time maps (panel b) show the supergranular patterns most 
clearly. 

FigureJTJi shows a 4-hour averaged intensity image in the 
Ca II H line of the north polar region. It is well-known that 
supergranulation is seen as a bright network structure in 
the chromosphere because magnetic field is concentrated 
at the supergranulation boundaries. In this chromospheric 
image the network brightenings, though faint, are seen, 
and their distribution is consistent with the supergranu- 
lar boundaries detected in the outward-inward travel-time 
maps (panels a to c). 

In Figure [TJ the outward-inward travel-time difference 
maps of the south polar region (panel e), a region near the 
east limb (panel f), and a region around the disk center 
(panel g) are shown. Figure [2] shows the nine outward- 
inward travel-time difference maps of the south pole ob- 
tained for each of the nine observing runs in March of 2010. 
In these maps, the supergranular patterns in both polar re- 
gions and the east limb region are seen as clearly as in the 
travel time map of the disk center. In the high-latitude 
region, however, the cells seem to align in the direction 
roughly from north to south, with some tilt. The 'align- 
ment' patterns exist in both of the north and south polar 
regions; Figure [5] suggests that the 'alignment' is probably 
a characteristic structure in the polar regions. There also 
seem to be areas with weak diverging flows between the 
aligned cells. The alignment seems to be coherent for a 
distance of 2-3 cells. 

4.2. Correlation of the travel-time difference maps 

For a quantitative characterization of the alignment, we 
calculated two-dimensional correlations of the travel-time 
difference maps. Figure [3] shows the correlation functions 
of the travel-time difference maps for the north polar re- 
gion and the lower-latitude regions as well. The correlation 
function of the north polar datasets (panel d) does indi- 
cate a striking alignment in the map, i.e., the correlation 
is the strongest in the northeast-southwest direction. On 
the other hand, for the disk center and the east limb region 
maps (panels e and f) the correlation shows no particular 
preference in direction; namely, it shows a nearly random 
distribution. 

Actually, the strength of the 'alignment' differs from 
dataset to dataset. The alignment, however, does exist in 
the averaged correlation functions of the travel-time dif- 
ference maps, calculated for different latitudes (Figure 2]). 
The alignment is most notable at the highest latitude, as 
shown in the bottom panel. The alignment is almost in 
the north-south direction, but not precisely. We estimated 
the mean direction and its variance for the nine datasets. 
We measured the peak position angle along a 20-Mm cir- 
cle centered at the origin of the correlation maps, and 
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found that the mean and the variance alignment angles 
are 90.1 ± 37.0 degrees at 72° in latitude, 102.5 ± 23.2 de- 
grees at 76°, and 92.1±11.5 degrees at 83°. Note that here 
the direction of degree is to the west, and 90 degrees is 
to the north. This analysis indicates: 1) the alignment is, 
on average, in the nearly north-south direction, and 2) the 
variance of the alignment angle is smaller in the higher lat- 
itude regions. On a short time scale, the alignment might 
be affected by t he polar dynamics, such as varying differ- 
ential rotation ()Ye &: Livingston 1998D . meridional flow, 
and a polar vortex (e.g., Gilman 19791 ). 

The correlation functions of the travel-time maps also 
show us variation of the cell size with latitude. We have 
found that the east-west size of the supergranular cells 
in the high-latitude regions tends to be smaller than in 
the low-latitude regions. We define the cell radius (a 
half of the full size) as a correlation length: the dis- 
tance in x (east-west) direction from the origin to the 
point where the correlation is down to zero in the cor- 
relation map. The cell correlation radii are 10-12 Mm 
in the north/south polar regions, and 13-14 Mm in the 
east limb region and around the disk center. Moreover, 
from the average of the correlation function over the 
nine south polar datasets, we obtained the cell radii as 
12.0 ± 0.4 Mm around 72° in latitude, 11.8 ± 0.2 Mm 
around 76°, and 11.2 ± 0.2 Mm around 83°. The er- 
rors are estimated from the variance of the nine datasets. 
This implies that the cell size depends on latitude. This 
type of latitudinal dependence is consistent with the re- 
sults obtained fr o m f-m ode helioseismology analysis by 
Hirzbcrg er et al.l (|2008f) and from Ca II K line obser- 
vatio n s of the chromospheric network dBrune fc Wohll 
Il982t iMunzer et al.q989)Munzer. Hanslmeierl ISchroterT h 
although these observations were for mid- or lower-latitude 
regions. Note that by this method we could not define the 
cell size in the north-south direction because the correla- 
tion function does not drop to zero in this direction (see 
Figure S]). This might be related to the fact that regions 
with weak diverging flows are located between the aligned 
cells in the north-south direction as can be seen in the 
travel-time maps (Figures Q] and [2]) , but this suggestion 
needs further investigation. 

We verified that the alignment is not due to some kind 
of artifacts caused by the limb observation, because the 
alignment is not found in the control study of the east 
limb region. In general, any kind of limb-observation ef- 
fect would not affect positions of the cells in the travel-time 
maps, thereby producing a spurious alignment. Also, it is 
unlikely that projection effects affected our estimates of 
the cell size in the cast-west direction. 



4.3. Temporal evolution 

The upper panel of Figure [5] shows the combined 
outward-inward travel-time difference map of the nine 
south polar datasets of March of 2010. This is a view 
from above the South Pole of the Sun. The first observa- 
tion of March 3 is the rightmost one, and the following ob- 
servations are located sequentially counter-clockwise. We 
calculated the longitudinal difference between t he observ- 
ing d ates by using the differential rotation rate ()Snodgrassl 
1984) at the central latitude of each observation. Note that 
the supergranulation pattern in Figure [5] looks fainter than 



in the individual maps (Figure [5]) because of overlapping 
of the different time intervals. 

The lower panels show the last two datasets, the tem- 
poral separation between which was only 15.5 hours, and 
their fields of view significantly overlapped. Many fea- 
tures in the overlapping region look similar; a clear exam- 
ple is along the line at angle of 240 degrees. Obviously, 
the cells and their alignment pattern survive at least 15.5 
hours, which is consistent with the known fact that the 
supergranular cell lifetime is about a day. For a detailed 
study of the temporal evolution we need more consecutive 
datasets. 

5. DISCUSSIONS 

We have found evidence indicating that: 1) supergranu- 
lar cells tend to align roughly in the north-south direction 
(with some tilt) in the polar regions; the alignment is co- 
herent over the scale of a few cells, 2) there may be areas 
with weak diverging flows between the aligned cells, and 
3) the horizontal (east-west) size of the cells decreases to- 
w ards higher l atitud es. 

iLisle et ail (|2004l) found observational evidence of the 
alignment of supergranular cells in low latitude regions 
below 60 degrees, and suggested that the alignment is con- 
trolled by giant cells, because converging flows at the giant 
cell boundaries are expected to align parallel to the rota- 
tion axis at low latitudes, as predicted by numerica l sim- 
ulations (|Brun fc Toomrell2002t iMiesch et al.ll2006[ ). The 
alignment reported bv lLisle et al.l l 2004 ) was found during 
8-day observations, while our observing runs were only 16 
hours long. Since the lifetime of a supergranular cell is 
about a day, the alignment found by them is somewhat 
statistical, i.e., it does not mean that supergranules are 
aligned at a given time, as was found by our study. Yet, 
in both cases the alignment may be a result of giant-cell 
flows organizing smaller-scale flow fields. What we essen- 
tially need for further study is to determine the lifetime 
of the organized structures, and to compare with the ex- 
pected giant-cell lifetime (about a month), as well as to 
study the relationship with the global magnetic structure 
of the Sun. Also, it should be interesting to compare the 
polar supergranulation pattern with the global supergran- 
ulation pattern by using the latest helioseismology instru- 
ment, Helioseismic and Magnetic Imager (HMI) onboard 
the Solar Dynamics Observatory. 

Perhaps, the fact that the cells in the polar regions 
seem to have areas with weak diverging flows between the 
aligned cells can be explained by the convective rolls. If 
there are elongated cells in the north-south direction with 
counter-rotating convective motions, the outward-inward 
travel-time difference m aps would s how these patterns. 
However, the theory fe.g. lBusselll970t ) predicts the convec- 
tive roll patterns in low-latitude regions of the convection 
zone below the surface, which remain to be detected. 
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Fig. 1. — Outward-inward travel-time difference maps of the north polar region on 2009 September 25, for travel distances: (a) 8.4 Mm, 
(b) 14.4 Mm, and (c) 21.1 Mm. Panel (d) shows a 4- hour averaged intensity image in the Ca II H line of the region. Panels (e), (f), and (g) 
show the outward-inward travel-time difference maps of the south polar region on 2009 March 7, near the east limb on 2009 September 27, 
and around the disk center on 2009 December 11, for the travel distances of 14.4 Mm. The gray scale used for the travel-time difference maps 
covers the range from —0.5 minutes (black) to +0.5 minutes (white). 
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Fig. 2. — Outward-inward travel-time difference maps of the south polar region observed for 9 periods in March of 2010. The observing 
dates are indicated in the lower right corner. Note that the March 16 map is noisier than others because of short observation period. The 
gray scale is same as in Figure [T] 
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Fig. 3. — Travel-time difference maps and two-dimensional correlation maps of the travel-time difference. Left column: Travel-time difference 
maps in (a) the north polar region, (b) the disk center region, and (c) the east limb region. These maps are Postel's projected maps, and the 
sizes of the field-of-view are the same. Note that (a) and (c) is only a part of the full field of view. The gray scale is same as in Figure [T] 
Right column: Two-dimensional correlation maps of the travel-time difference maps shown in the left panels. Panels (d), (e), and (f) are for 
the north polar region, the disk center region, and the east limb region, respectively. The gray-scale range of the cross-correlation function is 
from -0.5 (black) to +0.5 (white). 
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Fig. 4. — Averages of the two-dimensional cross-correlation functions over all the nine south pole datasets for the central latitudes: 72° 
(top), 76° (middle), and 83° (bottom). The gray scale is same as in the right panels of Figure [3] 
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Fig. 5. — Outward-inward travel-time difference maps of the south polar region observed in March of 2010. The upper panel shows the 
combined polar view above the South Pole, obtained by overlapping the results of the nine observing runs. The lower panels are the last two 
consecutive travel-time difference maps. The concentric circles are the latitudinal lines spaced by five degrees, while the straight lines are the 
longitudinal lines spaced by 30 degrees. The gray scales used for the upper panel and the lower panels cover the range from —0.25 minutes 
to +0.25 minutes and the range from —0.5 minutes to +0.5 minutes, respectively. 



